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1. Introduction 

In this paper, we investigate the quantum electromagnetic field in the vacuum region 
outside a half-space filled with a uniform, non-dispersive dielectric. We compute the 
expectations of the squares of the electric and magnetic fields in this region. We have 
two motivations for this. 

First, the problem is natural in the study of quantum optics. Indeed, other workers 
have already investigated some aspects of this situation, for example, the effects of 
a nearby dielectric on atomic transition rates (see, e.g., Khosravi and Loudon 1991, 
1992). Here we compute the expectations of the squares of the electric and magnetic 
fields, thus providing predictions of the Casimir-Polder force on an (electrically or 
magnetically) polarizable atom near the dielectric. These predictions test the ultraviolet 
renormalization of the the theory at a deeper level than do the transition-rate ones. 

Our main motivation, however, comes from the hypothesized uses of negative energy 
densities to fuel exotic general-relativistic and thermodynamic effects. Serious workers 
have considered the possibility that negative energy densities might give rise to "worm 
holes," "warp drives" and "time machines." Such predictions depend on being able to 
generate persistent negative energy densities. At present, the only way that this might 
be achieved within reasonably well-understood physics is via Casimir-type effects. In 
the original Casimir (1948) effect, for example, the energy density due to the quantum 
electromagnetic field between two perfect parallel plane conductors is predicted to be 
negative. It should immediately be remarked that this negative energy density has never 
been directly observed. [| Still, it is this prediction which has generated an enormous 
amount of theoretical work, because the possible consequences are so spectacular. 

We wanted to know what would happen to the prediction of negative energy 
densities if the plates were no longer idealized as perfect conductors. A realistic 
treatment of this would require a theory of the quantum electromagnetic field in 
inhomogeneous absorptive and dispersive media at finite temperature. Such theories 
are only now under development (see, e.g., Matloob et al. 1995), so it seems wise to 
consider as a first step the case of a non-absorptive, non-dispersive medium at zero 
temperature. Thus we shall consider the case of a half-space filled with a material of 
(frequency-independent) dielectric constant e. While the case of a perfect conductor is 
formally the limit e f oo of this, an imperfect conductor is not well represented by such 
a model with e finite. So we shall not be able to make any positive predictions about 
the behavior of real conductors. 

Still, our results are strong enough to bear on the case of conductors. We shall find 
that, for dielectrics, finite-e effects cannot be neglected, especially in computations of 

§ Laboratory experiments measure the force between the plates, that is, th component T zz of the 
stress-energy (Sparnaay 1957, 1958; Lamoreaux 1997, Bordag et al. 1998). The energy density is T tt . 
These two operators do not commute. There is a connection between them, in that the long-time 
average of the force is minus the gradient of the energy, but present experiments seem far from being 
able to measure T tt . This operator may as a matter of principle not be directly observable; see Heifer 
1998. 
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the electromagnetic contribution to the energy density. In the case of two parallel half- 
spaces, these finite-e corrections do not alter the attractive nature of the Casimir force, 
but may contribute a positive, separation-independent energy density which dominates 
the negative, separation-dependent, Casimir energy density. This strongly suggests that 
only after a careful treatment of the physics of real conductors will we know whether 
the perfect-conductor idealization is adequate for computing the energy density in such 
cases. 

It is not easy to say accurately and briefly why a finite dielectric constant should 
modify the energy density to this degree, because the physics is non-local and depends 
on quantum interference. The presence of a polarizable medium in a region alters the 
field operators, by causing reflection and refraction of modes at the boundary. If the 
geometry is particularly simple (a plane interface) and the reflection sufficiently idealized 
(a perfect conductor), one has a great deal of cancellation. Small deviations from these 
idealizations can potentially lead to large effects. This is beacuse the energy density 
and the squares of the field strengths are defined by ultraviolet -divergent integrals (and 
must be renormalized) . 

To explain the situation more quantitativley, we first review some aspects of the 
Casimir effect, and then discuss the idealizations that have been made and how they 
might be expected to be modified in a more realistic treatment. 

Between two perfect parallel plane conductors, one finds that the renormalized 
energy density is given by 

I rr*\ \ / 7~~">2 T~) 2 \ / . . 

(Too)rcn = ^ E +B ' rcn = ~ 720/4 ' ^ ' 

where / is the distance between the plates. That this is independent of position can 
be shown on invariance grounds (and relies on the ideal, perfect-conductor boundary 
conditions). However, the electric and magnetic fields are not position-independent; 
one finds 

~ 2 7r 2 Hc 7c 2 hc3 — 2 sin 2 (it z/ 1) 

{ )ren = " 720^ + 16Z 4 sin\nz/l) ( } 

/r2\ t 2hc n2fic3 ~ 2 sin2 (WO ,o\ 
{ )rca ~ 720Z 4 16Z 4 sm 4 (nz/l) [ ' 

at distance z from one plate. Near one plate, as z \ 0, we find the asymptotic forms 

<#>- ~ + ^ (4) 

In other words, the renormalized expectations of E 2 and B 2 both diverge near a perfectly 
conducting plate, but there is a perfect cancellation between the divergent terms, leaving 
only a finite result. 

Several commments on this are in order. First, the negative expectation of B 2 
occurs because it is a renormalized quantity, and means that the fluctuations of B 
are less than those of the Minkowski vacuum. Second, the divergences of {E 2 ) ren and 
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(B 2 )ren as z J, are not expected to be physical, but rather arise from the idealized 
boundary conditions used. A real conductor would not be well approximated by a perfect 
conductor within atomic distances, and probably not within its plasma wavelength. 
Thus the expressions are really only expected to be valid when one is sufficiently 

far from the conductor to neglect atomic structure and finite skin-depth. 

Still, one is led to ask what would happen if the antisymmetry between the divergent 
parts of (E 2 ) rcn and (B 2 ) rcn could be disturbed. Could one produce energy densities 
much greater in magnitude than the Casimir expression (D? A natural way to try to 
do this is to replace the perfect conductor by a dielectric, and this is what we have done 
here. Of course, our model is not expected to be accurate within atomic distances or 
even scales of the order of a skin-depth. Still, we shall be able to draw some interesting 
conclusions. 

We are able to compute (E 2 ) rcn , (E^) rcn , (B 2 ) ren and (B^) ien explicitly, as functions 
of the distance z from the dielectric boundary and of the dielectric susceptibility x- The 
expressions have the form 

r]hc/z 4 (6) 

where the 77s are transcendental functions of x- (See equations ^8], [34].) We find 



in particular that the energy density in the vacuum half-space has the form r] p hc/z 4 , 
where r\ p is a positive function of x- This means that the total energy per unit surface 
area of the electromagnetic field on the vacuum side, 

( Vp hc/z*)dz, (7) 







is divergent. This is unphysical and again can be ascribed to the oversimplification of 
our model, where all modes, of whatever frequency, are equally affected by the dielectric. 
In a more realistic model, the dielectric's atomic structure would be taken into account. 
This would mean that at small distances (of the order of the skin depth probably and 
at the atomic scale certainly) the energy density would not be given by r/phc/z 4 , but by 
some other, presumably finite, expression. Correspondingly, we ought really to think 
of our theory as an effective field theory valid only up to frequencies corresponding to 
wavelengths of order the skin depth or so. 

In the next section, we outline the technical details of the computations. In Section 
3, we summarize the asymptotic behaviors of the squares of the rjs for the squares of the 
fields, and present the graphs of these functions. Section 4 summarizes the behavior of 
the expectation of the stress tensor. Section 5 contains discussions of the significance of 
our results, and Section 6 recapitualiates the main conclusions. 



2. The Computation 

2.1. The Orthonormal Eigenmodes 

The case of a half-space uniformly filled with a dielectric has been studied earlier, and 
we shall use the orthonormal eigenmodes as given by Carniglia and Mandel (1971). 
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E E - kkT = (2 e )- 1 /2__ e ifc V ; ™(A;T x e z ) . (10) 



We shall take the z axis to be normal to the interface, with z increasing in the 
vacuum region. We take advantage of the translational symmetries in time and in the 
xt = (x, y) directions to resolve all modes by Fourier transforms in these variables, with 
Fourier transform variables uo and k^. These Fourier transform variables thus retain 
their senses on both sides of the interface. 

The dielectric constant is e = 1 + x- The wave number in the z-direction is k in 
the vacuum and k in the dielectric. Thus we have 

k 2 + k% = eu 2 for z < (dielectric) (8) 
k 2 + k\ = uj 2 for z > (vacuum). (9) 

In what follows k^ and e z are the unit vectors in the k^ and z-directions. In later 
sections, hats will indicate field operators, too, but no confusion should arise. 

We shall only need the modes on the vacuum side of the interface. The transverse 
electric component of the 'electric' field mode (where E is normal to the plane of 
incidence) incident from the left (k > 0), is 

2k 
k + k 

The transverse magnetic component of the 'electric' field mode incident from the left 
(k > 0) is 

2k / - \ 

E M ' kkT = (v / 2c)" 1 I e ^ e ^ _ kk \ . 

The transverse electric component of the 'electric' field mode incident from the right 
(k > 0) is 

E E kkT = V2~ l (V ifc + ^|e ifc ^ e ik —(k T x e z ) . (12) 

The transverse magnetic component of the 'electric' field mode incident from the right 
(k > 0) is 

E M kkT = (V2U)- 1 (e~ ikz (k T e z + kk^j + (^e 2 - kk^j ^ e ik ^ . (13) 

The transverse electric component of the 'magnetic' field mode incident from the left 
(k > 0) is 

2k 
k + k 

The transverse magnetic component of the 'magnetic' field mode incident from the left 
(k > 0) is 



B\ kT = (2eu 2 )^ 2 J^-e lkz e ik — (kk T - k T e z ) . (14) 



2k 

B M kkT = 2' 1 / 2 T ^-e ikz e ik —(k T x e z ) . (15) 
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The transverse electric component of the 'magnetic' field mode incident from the right 
(k > 0) is 

k — k 



B E kkT = -2- ll2 uo- 1 ( e~ ikz ( kk T + k T e z ) ^e ikz ( kk T - k T e 



(16) 



The transverse magnetic component of the 'magnetic' field mode incident from the right 
(k > 0) is 



B 



M 



fcfcx 



2 -l/2 e - ikz + 



ek — k 



jkz 1 a ik T -x T 



(fcr x e z ) 



k + ek 

The electric and magnetic field operators are thus given by 



;i7) 



?x "~ ia; * + h.c, 



W y fc>0 a=e,m v 



(18) 



A=E,M 



and 
S(x,t) 



(2 



n) Jk >° A=E,M 



t e~ iwt + h.c. 



+ 



(2 



^3 / d 3 A; J2 { a kk T Bkk T e ~ iuJt + h-c.) , 

/ J k>0 \_T? TV/T 



(19) 



A=E,M 

where the creation and annihilation operators satisfy the commutation relations 

4ir 3 h5x\i5{k-~k')5(k T -k T '), (20) 



a*' a^* 



and 



47r 3 M AA '5(A; - k')6(k T - fcr') 



(21) 



2.2. Computation of E 2 

We shall outline the computation of E^. Computations of the squares of the other field 
components follow the same pattern. 

We use a standard point-splitting in imaginary time, and set ir = t' — t. Then we 

have 



} (2vr)3 h 



3 ~k T 2 I 2k \ ( 2k 



d J k 



k -o 2w \ k + ek ) \ k + ek 

ek + k 



e i(k-k*)z e -UJT 



{2^fJ k>Q -2UJ 



e — + ^ e ikz e~ ikz 4 



e/c + A; 



We rewrite the integral over A; > as an integral over A; > (representing plane waves) 
plus an integral over < k < u^/x (representing evanescent waves): 



(£2) = JL f°° dco r dk(u 2 - k 2 ) [ 1 + ^4 cos2A;z 1 e" 



(2*0 2 7 7 



ek + k 
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o 



(2^) 2 

where we have used the relationships 



o.k(u + k)— — e e 



k 2 + e 2 K 2 



or 



+ for < cjy^x 
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(23) 
(24) 



uj 2 — k 2 for k > uj^Jx 

and have performed the simple polar angle integration. 

Changing variables k = u;£ in the first integral and k = ujy/xC, in the second integral 
we obtain 

4 



d^ 3 (i-a i+ 



et+Vx+e 



cos2ui^z e~ 



W? [ ^ /' d ^ 3(1 + Vxe)£0^e-^e- • (25) 



Integrating over uj gives 
he rl 



6— — + 6(1 -f )- 



+ 



126^(1 +x£ 2 )v / r r e 



(4^ 2 e 2 + r 2 Y 



(26) 



(l + (e 2 -l)e)(2^ + r) 4 _ 

Thus we have reduced the problem of finding the expectation value of the square of the 
^-component of the electric field to a one-dimensional integral. The integral can be 
evaluated using contour integration in the complex plane and by exploiting Cauchy's 
residue theorem. After integrating and extensive algebra we obtain the following 
formally divergent expression for the expectation value: 



lim 



he 



t^O V 7T 2 T 4 



+ 



he 



(2tt) 2 z 4 



16^( 2 ^ (6e2 - 3e3/2 " 2x) 



+6e(l-2e 2 + 2 X ) HV^+Vx) 

x-i) 



6e 2 ( 



2/^2 



In 



y/c + T- 1 



(Ve 



0{r) 



(27) 



Subtracting the (again divergent) vacuum (Minkowski space) expectation value 
(E 2 ) Minkowski = lim T _>o hc/ir 2 ^ and taking the limit as r — > gives the exact 
renormalized expectation value: 

+6e(l-2e 2 + 2 X ) HV~e + y/x) 
6e 2 fe 2 



{E% 



+ 



hcrr\ 



X 



In 



V7TT-1 
v^TT+ l 



(V7TT+ Ve~f 



(28) 
(29) 



say, where the coefficient rj^ is a function of x 
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The renormalized expectations of the squares of the other components can be 
calculated by the same techniques. They are given by: 

I f?2\ ^ C 
X-C'T/rei 



(-E^)ren 



and 
( b t)i 



(2tt)V 



16 ^ ( 2 Vx(6- 3^-2*)- 6(1 -2e*) 



ln(v^ + y/x) 



In 



Ve~+T+ 1 



(Ve + 1 + >/e) : 



(30) 
(31) 



(2vr) 2 2 



^(2^(12 -9v^-2 X ) 
-6(l-2 X )ln ^ +v ^)) 



(32) 
(33) 



(2tt)V 



1 



16x 3/2 



2 v ^(6 + 6e 2 -2 X -3( e + 2) v ^) 



+6(e-2e 3 + 2 X ) ln(^ + ^x) 

1 



+6e 2 ^e 



lln 



ye+T 
v^TT+ l 



(Ve + 1 + v 7 ^) 1 



(34) 



(35) 



These expressions are very complicated, and the characters of the functions rj(x) 
will be investigated in the next section. For the present, we remark that the sucessful 
renormalization provides a very strong check on the computations, since the term 
of order r~ 4 must cancel perfectly against the term from Minkowski space, and the 
remaining potential poles in r (of orders r~ 3 , r~ 2 and r _1 ) must vanish identically. 
Another check is provided by the vanishing of (T zz ) Ten , as will be discussed in Section 4. 



3. The Squares of the Fields 

In the previous section, we found the expectations of the squares of the fields explicitly. 
In each case the result had the form hcf]/z A , where z was the distance to the interface 
and r) was a complicated transcendental function of the susceptibility x- 111 this section, 
we present the graphs of the functions rj, as well as their limiting behaviors for x I 
and x T °° ■ 

The functions r] z ^ (scaled to have a common limiting value) are presented in 
figure 1. They are in each case monotonic and approach a constant value asymptotically, 
but the approach is extremely slow. 

For x I 0, we have 
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100 



Figure 1. The dependences of the squares of the fields on the susceptibility \ (the 
abscissa). The horizontal line at the top is the common asymptote. Below that, in 
descending order, are 7/:f , then (1/2)77^, then — (1/2)77^?, and finally — r]f. 



ftr 1 

< B «>'» = " (2^807* + (38) 

^>- = -^J^ X + °^- (39) 
It is good to observe that all the above expectation values tend to zero as x I (vacuum). 
For e ^> 1 (that is, x T we find 

^-"isfe-^i^ ^ (40) 

^~=w?-^w^ +0{m ' (41) 

( ^ = -16^ + (W^VJ +0(lM) ' ^ 

^» = -8^ + ^vl + 0(1/V?) ' <43> 
We note that in the limit x j 00, these quantities attain the values they would have 
in the case of the half-space outside a perfectly conducting plane (see e.g. Barton 
1990). This is in accord with the usual formal identification of perfect conductors with 
dielectrics of infinite susceptibility. However, the approach to this limit is rather slow. 
One needs x ~ 102 for rjf to be within 50% of its limiting value, and x ~ 14400 to be 
within 10%. 
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For completeness, we list the limiting behaviors of squares of the full fields: 
for x I 0, and 

, ^,oi 3 ftc f%c 231 _ . , . / , _ % 

< £ >™ = I6^- (2^807? + 0<1/£) <46) 

for x T °°- 

4. The Stress Tensor 

By symmetry considerations, the expectation of the renormalized stress tensor must be 
diagonal. The V component of the divergence constraint implies that {T zz ) rcn must be 
independent of z; however, as all components must be multiples of l/z 4 , this component 
must be zero. (The verification that one does get zero using our values of the squares of 
the fields provides another check on our calculation.) Since the tensor is trace-free and 
the 'xx' and L yy' components must be equal, there is only one algebraically independent 
component. We may take this to be the energy density p = (T tt ) ren , the other non-zero 
terms being (T xx ) Ten = (f yy ) rcn = p/(2c 2 ). 

The renormalized energy density is given by 

he 



(Too) i 



(2n) 2 z 



— (^(Oe 2 + 12 -A X - 3(6 + 3)Vi) 
+6 ( e _ 1 _ 2e 3 + 2(e + 1) X ) HV~e + Vx) 

(48) 



6eV-X-l),7v / f+T-l,/-^ rz\2 



say. 

The expectation of the energy density has the following limiting behaviours: 
for x ^ 1 we have 

™™ = (ik{r - 2 4) +o ^' < 50 > 

and for e ^> 1 we have 

The graph of the coefficient r] p is given in figure 2. It has a single maximum with 
value r] p pa .0034 at x ~ 13.65. The fall-off is again very slow, with i] p still at 50% its 
maximum value for x ~ 600, and still at 10% its maximum for x ~ 60 000. 
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100 



Figure 2. The dependence of the expectation of the energy density on the 
susceptibility. The coefficient r\ p is plotted as a function of the susceptibility \- There 
is a single niaxiniimum at \ ~ 13.65. 



5. Discussion 

We have found that the expectations of the renormalized squares of the components of 
the electric and magnetic fields at distance z from a non-dispersive dielectric all have the 
form r/hc/z 41 , where the coefficents rj depend on the susceptibility x- These coefficients 
vary linearly near x = 0, and tend monotonically but very slowly to asymptotic constant 
values as x T oo. For the electric field components, one has rj > 0; while for the magnetic 
field one has i] < 0. For the energy density, the corresponding coefficient i] p rises linearly 
near x = 0, attains a maximum at x ~ 13.65, and then falls very slowly to zero as x T °°- 
Previous authors have considered the quantum optical effects of a non-dispersive 
dielectric half-space on atomic transition rates (e.g. Khosravi & Loudon 1991, 1992). 
Those results are in some sense complementary to these: those test the field operators 
over a small range of frequencies; whereas the present ones depend on integrating 
over all frequencies. The present results depend crucially on the successful ultraviolet 
renormalization of the theory. 

5.1. Casimir-P older Forces 

An immediate consequence of these formulas is a prediction for the Casimir-Polder force 
on a polarizable atom near a dielectric. If the polarizability is cx(iv), then the induced 
dipole is f a{uo)E(uj) du>, and the potential energy is — (^f a(u)E(u) duo j ■ E. If in the 
regime in question we may neglect the frequencey-dependence of the polarizability, then, 
in the vacuum state of the electromagnetic field the potential energy is the vacuum state 
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of the electromagnetic field, this becomes 



- a{E 2 ) ren = -ahcr] E /z 4 , 



(52) 



where rj E = T] E + r] E . In principle, a parallel treatment applies to derive a 'magnetic 
Casimir-Polder' force depending on magnetic polarizability. 

The Casimir-Polder force on a polarizable atom near a conductor has been 
measured in recent years (Sukenik et al. 1993), although as yet there has been no 
measurement near a dielectric. 

5.2. Relation to Consitutive Energy 

One of our main motivations for studying this model was to uncover its limitations. We 
noted early on that a real physical dielectric cannot be well approximated as a uniform 
medium on arbitrarily small scales. This means that we cannot expect our model to 
accurately capture the physics of the field modes whose wavelengths are less than the 
atomic scale (certainly) or the skin depth (probably). In particular, our predictions 
must break down as one gets within a distance of this order of the interface. 

One vivid manifestation of this is the electromagnetic field energy of the vacuum 
half-space, per unit surface area. This surface energy density is 



which is divergent at the lower limit. (For an earlier investigation of a closely related 
effect, see Bordag and Lindig 1996.) It is not physically plausible that this energy 
density diverges; rather, the divergence reflects an improper model of physics near the 
interface. A more correct version would have the form 



where S is of the order of the skin depth. 

The energies we are considering here represent electromagnetic contributions to 
the constitutive energy of the medium. (In a more realistic treatment, it might not be 
meaningful to isolate one class of electromagnetic contributions from others, however.) 
It is plausible that these depend very much on the chemical physics of the material, and 
so the vagueness in the form ([54]) is apt. The numerical value of the contribution from 
the second term in ([54]) is quite modest for everyday materials. Taking the rather small 
value 5 « 10A and n p = .003, we find n p hc/(35 3 ) ~ 3 ■ 10~ 4 cal/cm 2 . 

5.3. The Fields in the Dielectric 

In this paper, we have only treated the electromagnetic field in the vacuum half-space 
outside the dielectric. This case would seem to be of more interest than the field within 
the dielectric. Still, it is natural to ask what would happen there. 




(53) 




(54) 
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In principle, techniques like ours should apply to compute the operators D = eE 
and B within the dielectric. The integrals involved are more difficult than those on the 
vacuum side, though. 

Aside from technical difficulties in evaluating the integrals, there is another point 
which must be considered in the dielectric region. In that region, the ultraviolet 
asymptotics of the two-point functions are different than in Minkowski space. (Because 
the two-point functions are singular on characteristics, and there is a different speed 
of light in the dielectric medium.) This means that the renormalization cannot be 
accomplished by subtracting the Minkowski-space vacuum quantities. One could 
presumably renormalize by subtracting the quantities associated to a uniform dielectric. 
If one does this, then the local energy density differs from that of Minkowski space by 
an infinite amount. 

We believe the resolution to this point is the same as that discussed in the previous 
subsection. One cannot accept the present model as an accurate picture of physics at 
all scales, and it is really only to be considered as an effective field theory, valid for 
frequencies below some cut-off. The difference in energy densities should be finite, with 
one contribution due to the effective field theory with a cut-off, and another due to the 
details of the chemical physics of the medium. 

5.4- Sign of the Energy Density 

One of our motivations for this work was to get a better understanding of the negative 
energy density occuring in the Casimir effect. 

The Casimir effect — corresponding to two parallel plane conductors — is formally 
the limit as e j oo of two parallel dielectric half-spaces. However, this identification only 
holds at the limit, and only in the sense that in this case the reflections and refractions 
of the field modes at the dielectric interfaces approach the perfect-conductor boundary 
conditions as e | oo. A real conductor, with finite conductivity, has a dielectric function 
which is significantly dispersive and absorptive, and cannot be modeled by a constant 
large real positive e. Thus our present model cannot make any positive quantitative 
predictions about the Casimir effect for conductors. 

However, we shall show that at least for dielectrics, the effects of a finite e cannot 
be ignored, and that in realistic situations it seems most likely that the total expected 
energy density, including separation-independent contributions from the half-spaces, is 
positive. This suggests strongly that we must investigate the real physics of conductors 
before we can conclude that the total expected energy density between the plates is 
negative. 

Consider two non-dispersive half-space dielectrics, of the same susceptibility, 
parallel and separated by a distance /. Then the total energy (per unit cross-sectional 
area) will have the form 

E tot = E OQ + r] 2 hc/l 3 . (55) 
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The function rj 2 (x) nas been computed by Lifshitz (1956). Here is (twice) the energy 
of either dielectric in isolation. 

For the energy density, let z with —1/2 < z < 1/2 be a coordinate normal to 
the interface planes. Then near either interface one expects the energy density to be 
dominated by the physics of that interface, and so to be ~ n p hc/{z ± Z/2) 4 , where 
r] p (x) is the coefficient we computed previously, equation fl4*8|). We shall write these two 
contributions as p%(z) and p 2 (z). The total energy density will be 

Ptot(z) = Pl(z) + p 2 (z) + p Cas (z, I) , (56) 

where pcas( z ,l) must on dimensional grounds have the form f{z/l)/l i 1 and be less 
singular at the interfaces than pi, p 2 . Indeed, the integral f_y 2 f(z/l)dz = r] 2 l/(hc) 
must be finite. Thus / has at most mild singularities at the interfaces. 

The form of the function / is at present unknown. As a very rough approximation, 
we shall assume it is constant in z/l, that is, the 'Casimir' contribution to the expected 
energy density is uniformly distributed between the half-spaces. We may ask if this 
value dominates the contribution p\ + p 2 from the dielectrics, that is, if p to t is positive 
or negative. This question can be answered by comparing our results with those of 
Lifshitz (1956), who computed the force of attraction of the two dielectrics. We find 
numerically that the energy density p\ + p 2 dominates the average energy density unless 
X ^ 39 000. In other words, if the 'Casimir' contribution to the energy density is 
distributed uniformly, the total energy density would be everywhere positive between 
the half-spaces unless \ could be made to exceed ~ 39 000. 

A real dielectric exhibits absorption and dispersion; we can expect our model to 
be valid at distances of order z if the dielectric susceptibility is (nearly) a real positive 
constant for several orders of magnitude of frequency bracketing c/z. There seems to 
be nothing in the Kramers-Kronig relations preventing this from holding for the sorts 
of values of x discussed above. Still, the scales are extreme enough that one wonders 
whether such susceptibility functions are more mathematical curiosities than physical 
possibilities. In other words, unless remarkable materials exist, with x( w ) approximately 
a real positive constant > 39 000 for several orders of magnitude of u>, it seems unlikely 
the total expected energy density anywhere between the dielectric half-spaces will be 
negative. 

We should like to emphasize that while the contributions from the dielectrics 
probably make the energy density positive, their separation-independence ensures that 
they do not alter the usual predictions of the attractive force between the dielectrics. 
Indeed, we have used Lifshitz's results in our argument, and we accept his values for 
the force. 

6. Conclusion 

Negative energy densities were first discovered in quantum field theory with Casimir's 
prediction of an attractive force between two parallel perfect plane conductors. Since 
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then, there has been considerable speculation on what the physical consequences of these 
negative energy densities might be. 

The present model was introduced as a first step away from the idealization of 
boundary conditions induced by a dielectric of infinite susceptibility. It has been chosen 
for its relative mathematical simplicity, and it is unrealistic in that it neglects dispersion 
and absorption. Still, we find that the effects of finite-susceptibility contributions go as 
e -!/ 2 i ne anc i can De ver y significant: in the range that our model is likely to be valid, it 
seems that these contributions can make the total energy density between two dielectric 
half-spaces positive, while preserving the attractive force found by Lifshitz. 

One cannot, from our model, draw any definite conclusion about the behavior 
of the energy density between two real conducting plates. What sorts of separation- 
independent corrections there are, due to finite conductivity, are at present unknown. 
But it does seem clear that we will only be justified in having confidence in theoretical 
predictions of the energy density between two real conducting plates if we take into 
account finite-conductivity effects. 
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